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ABSTRACT 



I 

Q Context. The detection of a narrow SiO line emission toward the young shocks of the L1448-mm outflow has been interpreted as a signature 
^ ' of the magnetic precursor of C-shocks. In contrast with the very low SiO abundances (<10"'^) derived from the ambient gas, the narrow SiO 
C/2 , emission in the precursor component at almost ambient velocities reveals enhanced SiO abundances of ~10"". This enhancement has been 
, ^ , proposed to be produced by the sputtering of the grain mantles at the very first stages of C-shocks. However, modelling of the sputtering of 
grains has usually averaged the SiO abundances over the dissipation region of C-shocks, which cannot explain the recent observations. 
Aims. To model the evolution of the gas phase abundances of molecules like SiO, CH3OH and H2O, produced by the sputtering of the grain 
^ ' mantles and cores as the shock propagates through the ambient gas. We consider different initial gas densities and shock velocities. 

Methods. We propose a parametric model to describe the physical structure of C-shocks as a function of time. Using the known sputtering 
0^ yields for water mantles (with other minor constituents like silicon and CH3OH) and olivine cores by collisions with H2, He, C, O, Si, Fe and 
CO, we follow the evolution of the abundances of silicon, CH3OH and H2O ejected from grains along the evolution of the shock. 
, Results. The evolution of the abundances of the sputtered silicon, CH3OH and H2O shows that CO seems to be the most efficient sputtering 
^N) ■ agent in low velocity shocks. The velocity threshold for the sputtering of silicon from the grain mantles is appreciably reduced (by 5-IOkm s"') 
' by CO compared to other models. The sputtering by CO can generate SiO abundances of ~10"" at the early stages of low velocity shocks, 
00 , consistent with those observed in the magnetic precursor component of L1448-mm. Our model satisfactorily reproduce the progressive en- 
hancement of SiO, CH3OH and H2O observed in this outflow, suggesting that this enhancement may be due to the propagation of two shocks 
^ _ with \'j=30km s" andv,=60kms" coexisting within the same region. 
• y—i . Conclusions. Our simple model can be used to estimate the time dependent evolution of the abundances of molecular shock tracers like SiO, 
■ CH3OH, H2O or NH3 in very young molecular outflows. 

■ - - Key words. ISM: clouds - physical processes: shock waves - ISM: jets and outflows - ISM: dust, extinction 



1. Introduction 

In young molecular outflows, it is expected that changes in tlie 
molecular emission could be observed due to the propagation 
of shocks into the ambient material. So far, the L1448-mm out- 
flow is the only object where time variability of the SiO emis- 
sion in the high velocity jet has been detected, indicating the 
presence of very young shocks (Girart & Acord 200J1). 

It is well known that silicon is heavily depleted onto the 
grain mantles and grain cores in the quiescent gas of molecu- 
lar dark clouds like TMC-1, L183 and L1448 (SiO abundance 
of <10 Ziurys, Friberg & Irvine 1989; Martin-Pintado, 
Bachiller & Fuente 1992; Requena-Torres et al. 2007). The de- 
tection of large SiO abundances in regions with outflow activity 
is therefore a clear indicator of the destruction of dust grains by 
the interaction of magnetohydrodynamic (MHD) shock waves 



(or C-shocks) with the ambient gas ([Martm-Fintado et al .1 19921: 
Flower et alii 1 9961; ICaselli et alJll997l) . 



The typical SiO abundances measured in the high ve- 
locity gas of young molecular outflows like in L1448-mm 
( Martin-Pintado et al.l 19921) are of >10 which implies an en- 
hancement by more than 6 orders of magnitude with respect to 
the SiO abundances measured in the quiescent gas. The recent 
detection of very narrow SiO emission at almost ambient ve- 
locities toward this outflow has been prop osed to be produced 



by th e magnetic precursor of C-shocks dJimenez-Serra et al 
2004 . The SiO abundance of -10 " for this narrow emis- 
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sion clearly contrasts with the large SiO enhancement found 
in the high velocity postshock gas, and with the much lower 
SiO abundance of the quiescent material. 

Toward the young sh ocks of the LI 448 -mm outflow, 
Jimenez-Serra et al.l ( l2005h have also reported an evolution- 
ary trend of the SiO and CH3OH abundances (methanol is 
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the most abundant molecule after H2O in the grain mantles; 
Tielens & Allamandolall987h to be enhanced from the ambient 
gas to the moderate velocity component, as if the grain mantles 
would have been progressively eroded by the recent interaction 
of low velocity shocks. 

Modelling of C-shocks that includes only the sputtering of 
grain cores, shows that an appreciable fraction of silicon ma- 
terial starts to be ej e cted from grains fo r Vj >25-30kms~' 



teriai starts to oe ei i 
dFlower et alj Il996[ 



Casern et alj Il997t iMav et alj I2OO0I) 



Although these models predict SiO abundances consistent with 
those observed in the postshock gas (-1 0^^-10"^), the sputter- 
ing of SiO from the cores cannot reproduce the lower SiO abun- 
dances of ~10 " found in the narrow precursor component of 
L1448-mm. 

Calculations of the sputtering yield of silicon by heavy 
atoms like C, O, Si and Fe on SiOi and olivine (MgFeSi04) 
cores show that, despite the low relative abundances of these 
species with respect to H2 and He in dark clouds, these heavy 
particles can dominate the sputtering of gr ains at low shock ve- 
locities ( Field et al.ll997l;lMav et al.l200Cil) . Furthermore, abun- 



dant molecules like CO could also play an important role in the 
sputtering of dust grains since these species can s putter like 
atoms of equivalent mass for low impact velocities (iMay et al. 
l200ft) . Considering that silicon could be a minor constituent 
of the mantles, their sputtering by these heavy species in low 
velocity shocks could efficiently erode them generating the 
SiO abundances observed for the narrow SiO line emission in 
LI 448 -mm. Up to date, the evolution of the sputtering of grains 
has not been studied in detail. The questions of which species 
are the most efficient sputtering agents, and which time-scales 
are needed to eject most of the silicon material from grains, still 
remain uncertain. 

In this paper, we present a parametric model of C-shocks to 
describe in detail the time dependent evolution of the molec- 
ular abundances sputtered from grains in low and high ve- 
locity shocks. This approximation constitutes a powerful tool 
for interpreting the molecular abundances measured in young 
molecular outflows. In addition to H2 and He, heavy atoms 
and molecules have been also considered as sputtering agents. 
In Sec. 2, we present the approximations used to describe the 
steady state profile of the physical structure of C-shocks. In 
Sec. 3, we show the procedure used to determine the sputtering 
of the grain mantles and the grain cores. In Sec. 4, we present 
the results of the sputtering of silicon from grains for several 
initial gas densities and shock velocities. In Sec. 5 and 6, we 
compare the sputtered SiO, CH3OH and H2O abundances with 
those measured in the LI 448 -mm outflow. The conclusions are 
finally summarized in Sec. 7. 

2. The C-shock Structure in the Preshock Frame 

We consider a plane-parallel C-shock that propagates through 
the quiescent gas with velocity Vj. As a first approximation, 
we have assumed steady state profiles for the evolution of the 
physical parameters in the shock. The validity of this approx- 
imation, versus more recent time-dependent modelling of the 
physical structure of C-shocks, will be discussed in detail in 
Sec. 4.1. 



Table 1. Initial gas phase abundances of He, C, O, Si, Fe and 
CO. 



Element/ 


Abundance" 


Molecule 


[«(X)/n(H)] 


He 


0.1 


C 


7.1x10-' 





1.8x10^* 


Si 


8.0x10-' 


Fe 


3.0x10-' 


CO 


1.5x10-* 



From the low metal case for dense cloud chemistry 
jPrasad & Huntress! 1 19821 : iHerbst & Leun3 1 19891 : lOraedel et al] 
Il982l) . 



The initial H2 density and temperature of the ambient cloud 
are no and Tq, respectively. Since one of the aims of this work 
is to directly compare our results with observations toward the 
young LI 448 -mm outflow, it is convenient to consider that the 
velocities of the ion and neutral fluids, v, and v„, are in the 
frame co-moving with the preshock gas. These velocities are 
approximated by: 



(vj - vo) 



cosh [(z - Zq)/z„j] 



(1) 



where z is the spatial coordinate and the z„/z, ratio gov- 
erns the strength of the velocity decoupling between the ion 
and neutral fluids, zo corresponds to the distance at which these 
fluids start to decouple (see Sec. 4.1 for details on how to esti- 
mate these parameters). An additional velocity, vo, also needs 
to be considered in the equations for v„ and v, in order to avoid 
infinite compression of the far downstream gas (see Eq. [3] be- 
low). Vo depends on the shock parameters and is defined as the 
final downstream velocity of the ion and neutral fluids in the 
frame of the shock (in the preshock frame, this final velocity 
would be - vq; see Sec. 4.1). As shown in Appendix A, this 
velocity is tightly linked to the shock and Alfven velocities, 
and va, through shock jump conditions. 

The ion-neutral drift speed is =|v„ - v',|, and the neutral 
fluid flow time is calculated as (see Eqs. lA.7l and lB.10l l: 



r dz 

J - v„ 



(2) 



From the principle of mass conservation, the neutral den- 
sity, n„, is given by: 



V, - v„ 



(3) 



The temperature of the neutral fluid, T„, is approximated by a 
Planck-like function as: 



[flr (z - zo)]* 



exp[(z- zo)/zt] - 1 



(4) 



where br is an integer, and ut and zt are related to the max- 
imum value of T„ iT„^max) and the distance Zn,max at which !„ 
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Table 2. Fractional abundances of H2O, CH3OH and Si/SiO 
assumed for the icy mantles and the grain cores. 



Species 


Abundances Ix] 




Mantles 


Cores 


H2O 


7.25xlO-5<"' 




CH3OH 






Si/SiO 




~3.6xlO-'<''> 



From'Whittet & Dulev (1991). 
From Jimenez-Serra et al. ( 2005 
f*^* From lAnders & GrevessJd 19891 



and lSnow & wit^ \l99dt) . 



reaches its maximum value. The temperature of the ion fluid is 

calculated by using T,- - T„+^^^. 

The comparison of the model predictions with observations 
(Sees. 5 and 6) requires the consideration of the radial velocity 
of the preshock gas (ambient cloud gas) relative to the observer, 
Vci- This velocity, the radial velocity of the emission measured 
by the observer, v^sr, and the velocity of the neutral fluid as 
measured in the frame of the ambient medium, v„, are related 
by: 



vlsr = Vcl + v„ 



(5) 



In Appendix A, we also give the equations for v„, v, and n„ 
within the frame of the shock (see Eqs. lA.ll and lA.5l l that will 
be used in Sec. 4.1 to validate this parametric approximation. 

3. Sputtering of Grains 

In this section, we describe the sputtering of grains produced by 
collisions with H2 and He and other heavy atomic and molec- 
ular species such as C, O, Si, Fe and CO (see Appendix B for 
the full explanation of the method). Although we consider that 
most silicon is locked into the olivine grain cores, we assume 
that a small fraction of this element is also present within the 
icy water mantles (^,,,=1.4x10"*; see below). CH3OH has been 
also considered as another constituent of the icy mantles. 

3.1. Sputtering of the Grain Mantles 

To study the sputtering of th e grain mantles, we have followed 
the procedure described bv ICaselh et al. ( 1997 ). The sputter- 
ing rate per unit volume and grain (Eq. lB.ll in Appendix B) 
has been derived by averaging the sputtering yield at low ener- 
gies (Eq. lB.2l ) over a velocity-shifted Maxwellian distribution 
characterized by T„ and v^. The surface binding en ergy Uq of 
the water mantles is of 0.53 eV (ITielens et al.lll994l) . The pro- 



jectile masses nip are 2, 4, 12, 16, 28, 56 and 28amu for H2, 
He, C, O, Si, Fe and CO, respectively. The target mass M, 
is considered to be ISamu which corresponds to the molec- 
ular mass of H2O. The initial fractional abundances of He, 
C, O, Si, Fe and CO, relative to atomic hydrogen, are shown 
in Table[T] We assume that these abundances remain constant 
throughout the dissipation region of the shock. The volume 



density of grains, rig, is derived by considering a gas-to-dust 
mass ratio of -100, a constant grain radius of 0.1 jum and a 
density of the grain core material of 3.5 gcm"^ (most of the 



volume of a grain is filled by the silicate core; ICaselli et al 



19971). The total sputtering rate for H2O, CH3OH and siU- 
con are calculated with Eqs. ( IB.4I ). ( IB.5l l and ( IB.6I 1. The total 
volume densities of these species are finally estimated with 
Eqs. ( IB.7I 1. ( IB. 81 ) and ( IB.9l l. Since the amount of material con- 
tained within the grain mantles is limited, we assume that the 
maximum abundances of silicon and CH3OH ejected from the 
mantles are those of SiO and CH3OH measured in the low ve- 
locity gas of the L1448-mm outflo w (~10~^ and ~1 0~^, respec- 
tively; see Table|2] and Jimenez- Serra et alJ l2005l) . From the 
observations , and assuming a H2O a bundance of ~7. 25x10"^ 
(Tablel2l and IWhittet & Dulevlll99lb . we can derive the frac- 
tion of silicon, q^, and CH3OH, r„„ present within the wa- 
ter mantles as q,„ = x(SiO)/x(H2O)^1.4xl0-'^ and r„, = 
X(CH3OH)/x(H2O)-lAxl0'^. Although q',„ and r„, constitute 
free parameters, we have fixed their values for comparison pur- 
poses with the L1448-mm outflow (Sees. 5 and 6). 



3.2. Sputtering of the Grain Cores 



For the sputtering of the cores, we have used different ap- 
proaches to calculate the sputtering produced by collisions with 
H2, and by collisions with He, C, O, Si, Fe and CO. We as- 
sume that olivine (MgFeSi04) is the main form of silicates in 
the cores. In the case of H2, we calculate the angle-averaged 
sputtering yield as in Sec. 3.1 (see Eq. B.2 in the Appendix B), 
but considerin g a surface bi nding energy Uo-SJOeW for the 



silicate cores (iTielens et al.l 11994) . For the sputtering agents 



He, C, O, Si and Fe , we have used t he sputtering yields for 
olivine calculated bv lMav et al. I tOOd see Eq.lOi. Since CO 
has the same projectile mass as Si, we ass ume that the sput- 



tering yield of CO is similar to that of Si dField et al.l 11997 



May et al.l 120001) . The sputtering threshold energies £,/, used 



for the projectiles are 73 eV for H e, 48 eV for C, a nd 47 eV 
for O, Si, Fe and CO (Table 4 of iMavet al.l l2000). To take 
into account the projection effects in the production of silicon 
within the shock, we have included the factor l/cos^^ in the 
impact energy Ep of the colliding particle (Eq. |B.3l ). where 6 
is the incUnation angle of the outflow with respect to the line 
of sig ht. For the L1448-mm outflow (see Sees. 5 and 6), is 
-70° jGirart & Aeordl200lb . As for the mantles, the total sput- 
tering rate for silicon is determined by Eq. ( IB. 6b . and (Eq. |B.9b 
calculates the volume density of silicon ejected from the grain 
cores. In the case of H2, we also need to consider that the prob- 
ability for a silicon atom to be injected into the gas phase from 
an o livine molecule, as opposed to a Mg, Fe or O atom, qc, is 
0.2 dCaselli et al.lll997l) . However, for the rest of the colliding 
particles, we assume that qc-^ since this probability has been 
already t aken into account in the calculations of the sputtering 
yields of lMav et alj j2000 ^. PracticaUy aU silicon (-99.97%) is 
locked into the grain cores with an abundance of ~3.6xl O~^ 
(see Tablel2t lAnders & Grevessi 19891; ISnow & "Mt3ll996l) . 
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Table 3. Input parameters for the C-shock profiles shown in Figs.[T]and|2] 





n(H2) 


So 




z„ 


Zi 


Zt 


Zo 


a-j- 


A 


T 

n,max 


(kms-') 


(cm-3) 


(mG) 




(cm) 


(cm) 


(cm) 


(cm) 


(K'/*cm-') 


(cm) 


(K) 


40 


10-* 


0.14 


7x10-*^ 


7.0x10'^ 


1.0x10'^ 


1.8x10'^ 


8.0x10'^ 


9.0x10-'^ 


5x10'" 


2000 


35 


10« 


14 


7x10-'° 


6.0x10'^* 


2.0x10" 


1.3x10'^ 


0.0 


1.3x10^''' 


2x10'" 


1700 




2 111'" 4 1(1" 

diNLan^e (cm I 



I I I I I I I I I I I I I I I I I I I I I I I I I I P M I I I I I 





(I 10° 2 iU" A 10'° 6 10" 

(lisUmw (cm) 



10' 



! I I I I ,i; I I I I I r 1 I M I I I I I I I I I — I I I I I I t I;: 




distance (cm) 



Fig. 1. Comparison between our approxi mation of v„, v,, Av, T„ and T, (lower panels), and the MHD shock structure calculated 



as in 



Flower & Pineau des ForetsI ( |200 3l upper panels) for a shock with V5=40kms «o=10"*cm ^ and Bq = lOOjuG (Flower, 



private communication). Shock parameters are given in Tableland velocities are in the shock frame (see Appendix A for the 
definition of v„ and v,). 



4. Results 

4.1. Validation of the Approximation of tlie Pliysical 
Structure of the C-shocl< 

To validate the parametric approximation of the C-shock phys- 
ical structure of Sec. 2, in Fig.[T|we show the comparison be- 
tween the C-shock profile calculated as in the recent MHD 
models of iFlower & Pineau des ForetsI (l2003h . and the profiles 
of v,„ V, , Av, T„ and T, derived in the shock frame through our 
approach (see Appendix A for details). In Fig.|2l we directly 
compare our results with the MHD shock structure obtained by 
Kaufman & Neufeld (Il996h . The values of z„, zi, Zt, Zo and aj 
chosen to reproduce these shock profiles (see below for the esti- 
mation of the input parameters), and the values of the magnetic 
field, Bq, fractional ionization, Xe, shock length scale. A, and 



maximum temperature of the neutral fluid, Tn,max, are shown in 
Tabled 



From Fig.[T] it is clear that, although some differences do 
exist between the approximation and the MHD shock mod- 
elling at moderate preshock densities (Flower, private commu- 
nication), the general behaviour of v„, v,, Av, r„ and T,, qualita- 
tively mimics the physical structure of C-shocks. In particular, 
the velocity decoupling between the ion and neutral fluids in 
the magnetic precursor, and the initial delay in the switch on of 
the neutral heating at this stage, are well reproduced by our ap- 
proximation. Note that the agreement between the profiles of v„ 
and V, is excellent, giving a reliable prediction of Av which is a 
key parameter in the sputtering yield calculation (see Appendix 
B). In order to fit the delay of the heating of the neutrals at the 
magnetic precursor stage, we need to impose bT>6. 
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0,5 1 1^ 1.5 

Distance [10 cm] 

Fig. 2. Comparison between the profiles of v„, v, and T„ de- 
rived with our approximati on (lower panel) and the M HD 
shock structure obtained by iKaufman & Neufeld d 19961) for 
a shock with Vs=35kms"' and «o=10**kms"' (upper panel). 
Parameters are given in Tableland velocities are in the shock 
frame. 



^ 10 ^ 



6 




0.006 



Fig. 3. C-shock physical structure obtained with Eqs. (1), 
(3) and (4) for v,=40kms-', no^lO^cm-^ and To^lOK. 
Velocities are in the frame co-moving with the preshock gas. 
The magnetic precursor length is of Az~0. 0005-0.001 pc=1.5- 
3.0x10'"^ cm. 



In Table m we summarize the different values of z„, z,, zt 
and aj, chosen to reproduce the physical structure of a sam- 
ple of C-shocks with velocities of 10<V5<40kms ' and ini- 
tial H2 densities of 10'*<n()<10^ cm"^. z„ and Zi have been esti- 
mated by considering that v„ is 0.999 Vj at A, and by assum- 
ing a Znizi ratio of ~9/2. We note that slightly higher (fac- 
tor of 1.5) ZnIZi ratio s are re quire d to reproduce the results of 
Iplower & Pineau des ForetsI ( 2003b with a new treatment of the 
coupling between the neutral and the charged fluids (Flower, 
private communication). However, for consistency, we will 
hereafter use the Znlzi ratio of ~ 9/2, since it well reprod uces 
the results of 'Flower et al.' ('1996') and lOraine et al.1 (Il983l) . for 
which Dopita & Sutherland (2003) accordingly give an esti- 
mate of the shock length scale, A. In any case, the results ob- 
tained with both values of z„/z, do not differ by more than 15%. 

The parameters aj and zt have been derived by assum- 
ing that i>7-=6 and zo=Ocm. The magnetic field, Bq, and frac- 
tional ioniza t ion of the gas, Xe, have been calculated as in 
Draine et al.l (1983 ). and th e shock length scale. A, as in 
Dopita & Sutherland! (l2003h . The estimated Alfven velocity 
is of VA=2.18kms"', and vo typically ranges from 3.1 to 
4.7 km s ' for the cases considered in Table 4 (see Appendix A 
for the calculation of va and vq). We note that the shock length 
scales derived for every initial H2 density of TablelH are of the 
same order of magnitude as the ion - neutra l coupling lengths 
determined by iKaufman & Neufeld (1 19961 see Fig. 1 in this 
work). The maximum temperature of the neutral fluid, T n.max, 
has been estimated from the results of iDraine et al.l (119831) . 



In Fig.|3] we show a representative profile of a C-shock 
with V5=40kms"', no=10^ cm"-' and Tq^IOK. As expected in 
the frame co-moving with the preshock gas, the ion and neu- 
tral fluids are at rest at the beginning of the shock, and their 
final velocities in the f ar downstream gas are of ~ Vj - vo 



(see iDraine et al.l Il983h . The initial delay of the heating of 



the neutrals gives the magnetic precursor length, which is of 
Az~0.0005-0.001 pc~1.5-3.0xl0'5 cm (Fig.©. While the max- 
imum value of the temperature of the ions is correlated with 
the maximum value of (see Fig.[3]l, the neutrals show their 
maximum temperature at Vn~0.8 5 (derived by assum ing H2O 
cooHng and a^^l.S in Eq. 18 of 'Smith & Brandl| l990l). wh ich 
is consistent with the results of Kaufman & Neufeld see 
Fig.H). 



For higher preshock densities, Fig.|2] shows that the ion- 
neutral velocity decoupling in the magnetic precursor, is not 
as ac curately reproduced as in Fig. [U fo r the moderate density 
case (IFlower & Pineau des ForetsI 120031) . We should mention. 



however, that the MHD treatment of C-shocks for high den- 
sity regions, is still rather simplistic and therefore, uncertain 
(see IPilipp et al.' 1990; Pilipp & HartquistI 1 19941; iFalld 12003 : 
Chapman & Wardle ,2006) . Since the existing perpendicular 



shock models may not start having p roblems until the p reshock 



density exceeds 10^ cm ^ (see e.g. Pilip p et al. 1990l) . in the 



following we will restrict our study to the moderate density 
case (from 10"^ to 10^ cm"^). 



From the recent results of time-dependent shock modelling, 
one may consider that our assumption of steadiness for the C- 
shock could not be valid to describe the time evolution of the 
sputtering of grains. These models indeed show that a J-type 
component is a natural feature in the far downstream gas of the 
C-shock (near the piston) for time-scales of <10^-10^ yr, for 



which the steady st ate is finally attained (IChieze et al.lll998 



Lesaffre et al. 2004). However, as shown in Sec. 4.2, the evo- 



lutionary stages relevant to the main injection of the material 
contained in the icy mantles and in the grain cores are those of 
the magnetic precursor which, independently on the age of the 
shock, can be described by the steady st ate profile of C-shocks 
dChieze et al.lll998HLesaffre et al.ll2004 . 
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Table 4. Input parameters for a sample of C-shocks. 



v., 




n(H2) 


Bo' 






z„ 


Zi 










y e 

^ n,max 


(km s-') 


(kms-i) 


(cm-3) 


(mG) 






(cm) 


(cm) 


(cm) 


(K'/*cm 


-') 


(pc) 


(K) 


20 


3.8 


10"* 


140 


7x10- 


-8 


1.4x10'^ 


3.2x10'^ 


5.0x10'^ 


2.9x10- 


16 


0.024 


900 


40 


4.7 


10-* 


140 


7x10" 


-8 


2.8x10'* 


6.2x10'-"^ 


1.1x10'* 


1.5x10- 


16 


0.048 


2200 


10 


3.1 


10^ 


450 


2x10- 


-8 


7.7x10'" 


1.7x10'" 


2.0x10'" 


5.8x10- 


15 


0.0012 


300 


20 


3.8 


10^ 


450 


2x10" 


-8 


1.4x10'^ 


3.2x10'" 


5.0x10'" 


2.8x10- 


15 


0.0024 


800 


30 


4.3 


10= 


450 


2x10" 


-8 


2.1x10'^ 


4.7x10'" 


8.0x10'" 


2.0x10- 


15 


0.0036 
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Calculated with Eq. lA.4l (see Append ix A) and assuming v^t =2. 1 8 km s . 
Estimated using Eqs. (62) and (63) of lOraine et al] ( Il983l) . 
Calculated c onsidering that fcj-=6 and zo= cm. 

Derived as in lOopita & Sutherland 1 20031) and ass uming that na Jnu ~10-*. 

Taken from Figs. 8b and 9b of lOraine et al.1 h983l) for ?Jci=10" cm and no =10^ and 10* cm ^ respectively. 



4.2. Sputtered Silicon Abundances: Injection and 
Saturation Times. 

We now include the evolutionary profiles of v„, v,, T„, T, and 
n„ from Sec. 4.1 in the sputtering equations of the Appendix 
B to calculate the silicon abundances ejected from the mantles 
and from the cores. Fig.|4]shows the silicon abundances ejected 
from grains as a function of the flow time for several H2 gas 
densities and shock velocities. The abundances of sputtered sil- 
icon do not practically change with the initial density of the gas, 
which is clearly in agreement with the results of Caselli et al. 
(Il997l) . However, as expected from the strongly dependence of 
die sputtering rate on the maxi mum value of vv (see Eg. IB. H and 



Pineau des Forets et al. 1997), the silicon abundance is dras- 



tically enhanced by increasing shock velocities. From Fig.|4] 
we also note that the time-scales are progressively reduced by 
nearly a factor of 10 as we increase the H2 density from 10^ 
to 10^ and 10^ cm"^. This is consistent with the fact that the 
flow time, f, is inversely proportional to the density (the cool- 
ing time-scales rou ghly vary as w,'', where w, is proportional 
to the density; see lChieze et al.il 1 9981; iLesaffre et al.il2004 . 

In Fig.|5] we show the products of the sputtering of the man- 
tles and of the cores for an initial density of 10^ cm and for 
shock velocities of 10, 20, 30 and 40 km s ' . The sputtering of 
the grain mantles by collisions with H2 and He (in bold lines 
in Fig.[5ll exactly co rresponds to that previously calculated by 
Caselli et al.l (1 1997b . Although the fractional abundance of the 
heavy species is orders of magnitude smaller than that of H2 
and He (see Sec. 3.1), it is clear that the heavy atoms and CO 
sputter the mantles much more efficiently than H2 or He for low 
shock velocities (Fig.|5]l. 

The high efficiency of these heavy species as sputtering 
agents is also shown by the injection, ti„j, and saturation times, 
tsat, of Table|5] We define f,„j as the time for which the gas 
phase silicon abundance, relative to H2, exceeds lO"^'' (i.e. the 
lower limit of Fig.|5]l; and f ,,,, as the time for which the dif- 
ference (in the logarithmic scale) of the silicon abundance be- 
tween two consecutive time steps f,+i and f,- (see Appendix A) 
is llogiolximjji+i - logiol)(:(m)]i\<0.1. Table|5]also shows the 



injection and saturation times due to the contribution of all col- 
liding particles. We note that these times are a factor of ~100 
smaller than the typical dynamical ages of young molecular 
outflows like L1448-mm (~1000yr), but are roughly of the 
same order of magnitude as the time-scales derived for the 
young shocks found in this outflow (<90 yr; iGirart & Acord 



2001b . From Table|5] we find that the injection and saturation 



times of H2 and He are larger than those of the heavy atoms 
and of CO. 

For the heavy colliding particles, CO seems to dominate 
the sputtering of the icy water mantles in low velocity shocks. 
Although Fe initiates the grain sputtering (Fe has the smallest 
injection and saturation times; see Table|5]), its low fractional 
abundance prevents large enhancements of silicon by the im- 
pact with this element (<3xlO for v., <20kms '; Fig.|5ll. 
On the contrary, collisions with CO (whose injection and sat- 
uration times are very similar to those of Si but whose initial 
abundance is 4 orders of magnitude larger than that of Si) pro- 
duce the main injection of silicon from the mantles. The satu- 
ration times for CO are indeed very similar to those derived for 
all colliding particles at low shock velocities (see Table|5l). 

The abundance of silicon sputtered by collisions with CO 
for v.5<10kms ' is very low (~10-'^; Fig.|5l). However, at 
slightly higher shock velocities (vj =20kms '), this abundant 
molecule can eject from the mantles considerable amounts of 
this element (~10-^). Averaging this silicon abundance over 
the dissipation region (shock length scale of ~7xl0'^ cm; see 
Tab.lU, we estimate that the total column density of silicon in- 
jected into the gas phase in a 20kms-'-shock is of ~10'^ cm"^. 
While shock velocities of Vj>25 km s ' wer e required t o obtai n 
Si/SiO column densities of >10'^cm"2 in Mav et all (l2000h . 
we find that the inclusion of silicon as a minor constituent of 
the grain mantles reduces the sputtering threshold velocity by, 
at least, | ~5 km s"'| in our model with respect to previous re- 
sults. This velocity thresh old is even reduced by | ~10kms"'| 
compared to the results of lCaselli et al.l (Il997b . 

It is also interesting to note that CO can also generate sil- 
icon abundances of ~2xl0"'' at the very early stages of low 
velocity shocks (f <10yr; see cases with Vs=20 and 30kms"' 
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Fig. 4. Evolution of the abundance of elemental silicon ejected from grains by the impact with H2, He, C, O, Si, Fe and CO, for 
initial H2 densities of 10^ 10^ and lO*" cm and shock velocities of 20 and 40 km s . 

n(H2)=10^cm"^ 
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Fig. 5. Abundance of elemental silicon ejected from grains by the impact with H2, He, C, O, Si, Fe and CO, for shock velocities 
of 10, 20, 30 and 40kms"' and a H2 gas density of 10^ cm"^. For each shock velocity, we show the individual production of 
silicon from the mantles (upper panels), from the cores (middle panels), and the total production of silicon from grains (lower 
panels). 



in Fig.|5]l. As discussed in Sec. 5, these results could explain the 
detection of SiO abundances of ~10 " associated with the nar- 
row SiO emission observed in the young shocks of the L1448- 
mm outflow. 

For shocks with v, >30kms"\ He plays an important role 
in the sputtering of the mantles. Note that the saturation times 
for all colliding particles at these shock velocities, slightly de- 
viate from those of CO due to the increasing efficiency of He to 
erode the grain mantles (Table|5]and Fig.|5]l. Almost all silicon 
within the mantles (~8xlO^^) is injected into the gas phase by 
collisions with He for shock velocities of ~30kms 

For the sputtering of the grain cores, only collisions with 
O, Si, Fe and CO are efficient enough to destroy the cores. The 
injection of silicon into the gas phase from the cores occurs for 



shock velocities Vj^30km s \ which is consistent with the re- 
sults of lMav et all ( I2OOOI) . For = 40 km s ', only -3% of the 
total amount of silicon locked into the grain cores, is released 
into the gas phase (abundance of ~10"*; see Fig.|5]l. As for the 
mantles, and although Si has the same injection and saturation 
times as those of CO (both are assumed to have similar sput- 
tering properties for the cores; see Sec. 3.2), CO is the main 
sputtering agent of the grain cores since its initial fractional 
abundance clearly exceeds that of Si. 

5. Comparison with Observations: The SiO 
Abundances 

If we now assume that silico n is rapidly oxidized into SiO 
jPineau des Forets et alJll997l) or that SiO is directly released 
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Table 5. Injection and saturation times for the grain mantles and the grain cores for a medium with an initial H2 density of 

no = 10^ cmT^ . 
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NOTE .- The injection and saturation times are given in yr. 



from grains dMartm-Pintado et al we can directly com- 

pare our predictions of the silicon abundance ejected from 
grains by the sputtering, with the SiO abundances observed in 
very young bipolar outflows like in L1448-mm. 

In Fig.|6] we show the SiO abundances measured for the 
different velocity components detected in this outflow (ambi- 
ent gas, the shock precursor compo nent, the moderate veloc- 



ity gas and the high veloc ity gas; IJimenez-Serra et al.l 12005 



Martm-Pintado et al.lll992h as a function of the flow time. For 
the ambient and precursor components, we have assumed cen- 
tral radial velocities of visr =4.7 and 5.2 km s ', respectively 
dJimenez-Serra et al.ll2004.) . Subtracting the ambient cloud ve- 
locity of Vci =4.7 km s"' characteristic of the molecular emis- 
sion in L1448-mm from the central radial velocities of these 
components (see Eq.|5] in Sec. 2), we obtain flow velocities 
of v„ =0 and 0.5 km s ', which correspond to flow times of 
f=0 and 4.0 yrs for the ambient gas and the precursor com- 
ponent in a 30 km s '-shock (Eq.|2]of Sec. 2; see also Eq. IA.7l 
in Appendix A for the details on the computation of the flow 
times). For the moderate velocity gas, we have considered 
velocity intervals of 1 km s '-width between 6 and 18kms"' 
dJimenez-Serra et al.ll2005h . The SiO abundances for the am- 
bient, precursor and moderate velocity gas have been derived 
assuming opti cally thin emission and exci tation temperatures 
of -10-15 K dRequena-Torres et al.1 l2007l) . For flie high ve- 
locity gas, the SiO abundances o f ~10~^ and 2x10"^ ( Fig.[6] l 
have been taken from Tab. 6 in iMartm-Pintado et al.l (Il992h 
for the velocity ranges of -50km s~'<VLSfl<-40km s~' and 
-60km s"'<yLSR<-50km s"'. The flow time associated with 
each velocity interval is again estimated from their central ra- 
dial velocities, vur, after subtracting the ambient cloud veloc- 
ity of the L1448-mm outflow of Vd=4.7 km s"' (Eq.|5]l. 

Fig. |6] also shows the silicon abundances, as a function of 
time, predicted by our model for the shocks which best fit the 
observational SiO data (with V5=30kms"' and Vj=60kms"'). 
The shock parameters used to reproduce the physical structure 
of these shocks are shown in Table|6] 

From Fig.|6] we note that the sputtering produced by 
the propagation of a 30 kms '-shock perfectly matches the 
evolutionary trend of SiO to be enhanced from the ambi- 



ent to the moderate velo city gas observed in L1448-mm 
(iJimenez-Serra et al h005'). The progressive erosion of the icy 
mantles by the sputtering with CO, reproduces the SiO abun- 
dances observed in the ambient gas (<10"'^; filled square), 
in the precursor component (~10"; filled circle), and in the 



moderate velocity gas (from ~10" to ~10" , filled triangles 
in Fig.|6ll. This suggests that the puzzling narrow SiO line de- 
tected toward the young shocks of L1448-mm can be explained 
by the recent erosion of the grain mantles containing a small 
fraction of Si/SiO, at the early stages of low velocity shocks. 

To fit the SiO abundances measured in the high velocity 
gas, a velocity shock with =60 kms"' is needed to sput- 
ter ~9% of the silicon contained within the olivine cores and 
increase the predicted silicon abundance up to a few 10"* 
(filled stars in Fig.|6]l. This shock velocity is clearly in ex- 
ces s of the critical velocities of C-shocks {vc r it ~40-50km s~'; 
jDraineetal."l983'; 'Smith & Braridlll990l). iLe Bourlot et al 



(I2002h and Cabrit et al, (2004) have recently shown that the 
actual shock velocity limit can be increased to v„/r~100km s"' 
for moderate densities and high magnetic fields. However, we 
cannot rule out the possibility that a J-type component would 
be the responsible for the large SiO abundances observed in the 
high velocity gas of L1448-mm. 

Given the fact that the L1448-mm outflow shows variability 
in its high velocity SiO emission, an alternative scenario would 
involve the presence of two different shocks at two different 
evolutionary stages that woul d coexist within the single -dish 
beam of the SiO observations dJimenez-Serra et al ] |2005h . 



6. Sputtered CH3OH and H2O Abundances. 

In addition to SiO, CH3OH and H2O are also expected to 
be largely enhanced in outflow r egions (see e.g. iDraine et al. 
19831 iKaufman & Neufeidlll996l) . In Fig. |7] (upper and mid- 
dle panels), we show the predicted abundances of CH3OH and 
H2O as a function of the flow time, for a sample of shocks 
with v,=10, 20, 30 and 40km s"' (see TableHli. The impact 
with CO injects CH3OH and H2O abundances as large as ~10"^ 
and ~10"^, respectively, in low velocity shocks (see cases with 
Vj=20kms"'; Fig.|7]i. These abundances are even enhanced to 
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Table 6. Parameters of the C-shock models that best fit the SiO and CH3OH observational data shown in Figs.|6]and|7J 
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Fig. 6. Predicted SiO abundances within a 30kms and a 
60 km s"' -shock as a function of the flow time (see shock pa- 
rameters in Table|6l). Observational SiO abundances derived 
in the ambient gas (filled square), the precursor component 
(filled circle), the moderate velocity gas (filled triangles) and 
the high vel ocity regime (filled stars) found in the LI 448 



mm outflow jMartm-Pintado et al.l Il992t iJimenez-Serra et al 



2005h . are also shown. The black arrow indicates an upper limit 
to the SiO abundance. The observed flow times have been de- 
rived from Eq.|2]of Sec. 2 (see text and Appendix A for details). 



As for SiO, the abundance of sputtered CH3OH is in agree- 
ment with that derived in the precursor component within one 
order of magnitude, and with those measured in the moderate 
velocity gas within a factor of <5 (see Fig.|7]l- We would like 
to stress that we have only argued the total amount of mate- 
rial in the grain mantles. This implies that the progressive en- 
hancement of the SiO and CH3OH abundan ces observed for ve- 
locities of <20kms ' toward L1448-mm dJimenez-Serra et al 



20051), is naturally explained by the presence of a single shock 
with Vj=30kms '. In the case of ortho-H20, however, the 
abundances of this molecule in the moderate velocity gas dif- 
fer by more than a factor of 10 from those predicted by our 
model. This could be due either to the assumption of op- 
tically thin emission for the ortho-H20 lio— >l()i transition 
(which could have und erestimated the ortho-H20 abundances; 
Benedettini et alJl2002 ), or to a beam dilution effect. Note that 
the SWAS beam is of -240", i.e., 8 times larger flian flie 30 m 
beam of the SiO 7=2— >1 observations (iJimenez-Serra et al 



2005h . This leads to a corrected ortho-H20 abundance of some 
which is consistent with our model predictions of Fig.|7l 



-10- 



up to -10-* for CH3OH, and to -10-"* for H2O, in shocks with 
only Vj=30km s-'. 

We can now compare the predicted abundances of CH3OH 
and H2O with those observed in the young shocks of the 
L1448-mm outflow. Fig.|7] (lower panel) shows the SiO and 
CH3OH abundances observed in the ambient gas (square), 
the precursor component (circles) and the moderate ve- 
locity gas (triangles) of this outflow (IJimenez-Serra et al 



20051) . as a function of the flow time. In this figure, we 
also show the abundances of ortho-H20 derived from the 
line profile of the lio ^lpi transition measured by SWAS 
dBenedettini et al ] l200i . Si nce the velocity re solution of these 
observations ( ~ 1 km s~' : iBenedettini et a D I2002) is lower 
than those of the SiO and CH3OH spectra (~0.14kms-'; see 



Jimenez-Serra et al.l 120051) . we have only estimated the abun- 
dances of ortho-H20 for the moderate velocity regime (dark 
grey triangles in Fig.|7]i. As in Sec. 5, the flow times associ- 
ated with these abundances have been inferred subtracting the 
ambient cloud velocity in L1448-mm of Vri=4.7kms-' from 
the central velocities, visr, of the observed line profiles of 
these molecules (see Eq.|5]l. The SiO, CH3OH and H2O abun- 
da nces have been derived as suming optically thin emission as 
in Jimenez-Serra et akl (1200 5) for SiO and CH3OH, and as in 



Neufeld et alJ (l2000h for H2O (Eq. 1 in this work). The abun- 
dance of SiO (black line), CH3OH (light grey line) and H2O 
(dark grey line) in gas phase generated by the sputtering of the 
grain mantles in the 30 km s-' -shock of Sec. 5 (see Table|6]l are 
also shown in Fig.]?] 



7. Conclusions 

In this work, we have presented a parametric model that mimics 
the steady state profile of the physical parameters of C-shocks. 
The simplicity of this model has allowed, for the first time, the 
detailed analysis of the time evolution of the sputtering of the 
grain mantles and the grain cores in regions with recent outflow 
activity. Although this approximation does not include detailed 
MHD modelling, we have shown that it can be used as an effi- 
cient tool to interpret the time dependent evolution of the abun- 
dances of typical shock tracers like SiO, CH3OH, H2O or NH3 
in young molecular outflows, where transient phenomena are 
expected to play an important role. The assumption of steadi- 
ness of the shock can be applied for these young objects, since 
only the early stages of the shock evolution (characterized ex- 
clusively by the interaction of the magnetic precursor) are rel- 
evant for the process of the sputtering of dust grains. 

To calculate the sputtering of the grain mantles and the 
grain cores, we have assumed that silicon and methanol are mi- 
nor constituents of the water mantles, and that olivine is the 
main form of silicates within the grain cores. In spite of the 
low fractional abundance of heavy atoms (C, O, Si and Fe) 
and molecules (CO) relative to H2 or He in molecular dark 
clouds, these species have been also considered as sputtering 
agents. The relatively high abundance of CO with respect to the 
rest of heavy colliding particles, and its large sputtering yield, 
makes the sputtering by CO very efficient. Collisions with this 
molecule can eject a considerably large fraction of silicon in 
the mantles for shocks with only Vj ~20kms-'. This implies 
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Fig. 7. Upper and middle panels: Predictions of the sputtered 
abundances of CH3OH and H2O as a function of the flow time, 
for a sample of C-shocks with Vj=10, 20, 30 and 40kms ' 
and «o=10^ cm"^. Lower panel: SiO, CH3OH and H2O abun- 
dances ejected from the grain mantles for the 30km s '-shock 
of Table|6]as a function of time. We also show the SiO, CH3OH 
and H2O abundances measured in the ambient gas (square), 
the precursor component (circles) and th e moderate velocity 
gas (triangles) of the L1448-r nm outflow djimenez-Serra et al 



2005t iBenedettini et alj|2002h . Vertical arrows indicate upper 
limits to the SiO abundance. 



a reduction by | ~5-10kms"'| of the threshold velocity of the 
sputtering with respect to other models. 

By comparing the evolution of the abundances of SiO, 
CH3OH and H2O predicted by our model, with the abun- 
dances derived toward the different velocity regimes found 
in the young shocks of the L1448-mm outflow, we find that 
two different shocks (with y 5=30kms ' and V5=60kms"') are 
needed to reproduce the measured abundances. The progres- 
sive enhancement of SiO and CH3OH observed from the ambi- 
ent gas up to moderate velocities, is consistent with the mantle 
erosion produced by a single 30 km s '-shock. We find that the 
SiO abundance of ~10 " associated with the very narrow SiO 
emission detected in this outflow, can be explained as an early 
product of the sputtering of the mantles by CO in low velocity 
shocks if the Si/SiO abundance in the grain mantles is ~10 
with respect to water The disagreement between the predicted 
and the derived ortho-H20 abundances is probably due either 
to the assumption of optically thin emission for the ortho-H20 
1 10— *loi line observed by SWAS, or to a beam dilution effect. 



The approximation presented in this work, not only will 
furnish an input for more comprehensive MHD models of the 
shock structure in young molecular outflows, but will allow to 
perform direct comparisons with the molecular line profiles ob- 
served toward these regions. These comparisons will be pre- 
sented in a future paper. 
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Appendix A: The C-shock Structure in the Shocit 
Frame 

The steady-state velocity profiles of the ion and neutral fluids, 
V* and V*, within the frame of the shock, are approximated as: 



V'o + 



(Vj - vo) 



cosh \{z - zo)/z«,i] 



(A.l) 



where z is the spatial coordinate and is the velocity of the 
shock, zo, Z\ and z„ are input parameters that govern the veloc- 
ity decoupling between the ions and the neutrals (see Sees. 2 
and 4.1), and vo is required to avoid infinite compression of the 
far downstream gas (see Eg. |A. 51 and Sec. 2). This speed only 
depends on the speed of the shock, v,, and the Alfven speed, va 
(see Eq. lA.4l below). From Eq. lA.ll it is clear that the ion and 
neutral fluids initially move with velocities Vj and are progres- 
sively decelerated to vq (see Figs.[T]and|2|i. 

Assuming that the thermal pressure of the fluid is negligible 
compared to the magnetic and the dynamic pressures, we can 
estimate the magnitude of vq in the downstream gas from: 



(A.2) 



where po and p are the preshock and postshock gas mass 
densities, and Bq and B are the magnitude of the magnetic 
field in the preshock and postshock regimes, respectively. 
Considering that the compression of the magnetic field is given 
hy B - Boivjvo) and that the Alfven speed, va, is defined as 
vl = Bl/iAnp), Eq.lA!2lbecomes: 



P/Po 



1-1^ 

Po 



(A.3) 



Jimenez-Serra et al.: Evolution of the Sputtering of Grains 



11 



From the principle of mass conservation, the fluid density 
in the postshock gas, p, is derived as p = Po(Vs/vo). If we now 
substitute the ratio p/po in Eq. |A.3| we finally obtain: 



Vq - VqVj 



(A.4) 



from which vq can be calculated. 

The particle density of the neutral fluid, n„, is given by; 



(A.5) 



and the time associated with the neutral fluid or flow time, 
t, in the frame of the shock is calculated as: 



I V 



In each plane-parallel slab of material / + 1 within the 
shock, the flow time, f,+i , is calculated by using the trapezoidal 
method: 



ti+l = ti + (Zi+l - Zi) 



fiZi) + f{ZM) 



(A.7) 



where z,- and Zi+\ are the spatial coordinates for the slabs of 
gas i and ; + 1, and /(z,) is defined as (see Eq. |A.6l ): 



Y(E) >g^2Y(E,6 = 0) ( lDrainelll995b . where the normal- 
incidence yield Y(E, - 0) is calculated as (Eq. 3 1 ; 
Draine & Salpeteij|l979h : 



Y{E, = 0) = A 



1 + (e/30)' 



4/3' 



E > £() 



(B.2) 



A is a constant (A x 8.3 x 10"^), and s and so are cal- 
culated as e = TjEplUd and sq - max[l,4r]]- Uq is the bind- 
ing energy (per atom or molecule) and 77 is derived by doing 
T] = A^nipMiinip + MtY^, where M, is the target mass and 
^ is an efficiency f actor which varies from 0.8 (for ices) to 1 
(for atomic solids; IPraine & Salpeterl[l979l) . x,h, which is re- 
lated to the threshold impact energy £,/„ is finally calculated as 



\ nkT„ I 



For the grain cores, we h ave used the sputtering yield cal- 
culated bv iMav et al. (2000') for the impact of atomic species 
on olivine cores. The sputtering yield is derived from: 



(A.6) < Y{E) >fl= k, exp [-(HiEp - E,h)\ 



(B.3) 



where fe,, fi and E,h (the sputtering threshold energy) are 
taken from Table 4 in lMay et~al. (.2000.) . In this case, jc,/, is de- 
rived by doing x,h — {jj-)- 

In each collision between projectile and grain, only a small 
fraction of silicon, q„i, and CH3OH, r,„, will be ejected from 
the mantles [q'„,= 1.4x10"'* and r„,=1.4xl0"^; see Sec. 3.1]. 
Analogously, only a fraction of silicon, qc, will be released 
from the cores {qc-Q.2 for H2 and qc-l for the rest of col- 
liding particles; see Sec. 3.2]. If we assume a grain density Ug, 
the total sputtering rate for H2O, CH3OH and silicon is: 



fiZi) 



Kizi) 



(A.8) 



dn(H20) 



dt 



dn(H20) 



dt 



(B.4) 



The temperature of the ion and neutral fluids, T, and T„, are 
estimated as in Sec. 2. 



dniCHiOH) 



dt 



dniCHiOH) 



dt 



gram 



(B.5) 



Appendix B: Sputtering of Grains. Silicon, Cl-laOI-l 
and H2O Fractional Abundances. 

The sputtering of grains has been calculated by considering dif- 
ferent sputtering yields for the mantles and for the cores. The 
sputtering rate per unit volume for a spherical target of radius 
a moving with drift velocity th rough a Maxwellian neu tral 
gas of temperature T„ is (Eq. 27 in lDraine & Salpete 



[dnim) 1 
~dr\ 



gram 



■ f dxx^ h\e 



dn(S i) 



dt 



— qm,c 



dn(S i) 



dt 



(B.6) 



where m and c denotes the sputtering rate for the mantles 
and the cores. 

By using the Euler's algorithm, we calculate the total vol- 
ume density of H2O, CH3OH and Si ejected from grains in each 
plane-parallel slab of material / within the shock as: 



(B.l) n(H20)M = n(H20)i + At 



] < Y(E) >e 



dn(H20) 



dt 



Jtot,i 



where rip and nip are the number density and mass of the 
projectile, respectively; and s and x are related to Vd, T„ and the 

projectile impact energy Ep through s-^ ™d Ep - x^kT„. 

The angle-averaged sputtering yield at low energies 
< Y{E) >g for the mantles can be approximated by < 



n{CH^OH)i^i = n{CHj,OH)i + At 



dniCHjOH) 



dt 



n(S i)i + At 



dn(S i) 



dt 



(B.l) 



(B.8) 



(B.9) 



tot,i 
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where Af = (?,+[ - f,). In this case, the flow time at the slab 
of material / + 1, f,+ i, is numerically calculated as in Appendix 
A (see Eq. |A.7| ). but using the function /(z,) defined as (see 
Sec. 2): 



f(Zi) 



1 



Vs - V„(Zi) 



(B.IO) 



The fractional abundance of silicon, H2O and CH3OH, is 
finally derived by doing x(m) - n(m)/n(H2) in each slab of 
material. 
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